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VOLUME 8, 1941 


SOME REMARKS ON THE MOVEMENT OF CHROMOSOMES 
DURING CELL DIVISION 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


The possibility of drawing conclusions about the nature of the 
forces acting upon the chromosomes during division from observa- 
tion of the rates of their movement in anaphase is pointed out. Some 
available data are discussed, and shown to agree quantitatively with 
the assumption that during anaphase the chromosomes are pulled 
apart by contracting elastic fibers. 


In a previous outline of a physicomathematical theory of mitosis 
(Rashevsky, 1938; chapter xiii) it has been suggested that the chromo- 
somes are kept in a stable position in the equatorial plate during meta- 
phase by repulsion diffusion forces from two “equivalent” diffusion 
centers at the poles of the spindle. During anaphase they become con- 
nected by elastic fibers with the poles. When the elastic forces exceed 
both the diffusion forces as well as the forces that may hold together 
the daughter chromosomes, these are pulled by contracting fibers 
towards the poles. 

Since the chromosomes move in a viscous medium and their ve- 
locities are small, the inertial terms in the equations of their motion 
may be neglected and we may set their velocity at any moment as 
proportional to the force acting at this moment and inversely pro- 
portional to the viscosity. Denoting by y the distance of a chromosome 
from the pole, the elastic force pulling it towards the pole will be 
approximately proportional to y, while the diffusion force, repelling 
it would be approximately proportional to 1/y?. The equation of mo- 
tion of a chromosome will be of the form: 


i (1) 
a and b being constants. Integrating equation (1) and denoting by 
y, the value of y for t=0, we find: 
ay? — 0 
ee ee te 9 
log eb 3a (2) 
For very small values of b, or very large values of y, equation 
(2) reduces to 
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log y = — at + Const. (3) 


In the presence of only elastic forces a straight line relationship would 
hold between log y and t. In the presence of repelling forces also, the 
more complex equation (2) will hold. 

Data that can be used to verify one of the two above equations 
have been supplied by H. N. Barber (1939), who gives relative dis- 
tances of chromosomes during anaphase for various cells. Unfortu- 
nately, not all of his data can be used for the present purpose. The 
three sets of curves for Tradescantia do obviously not satisfy the 
simple conditions considered here. In those cases the chromosomes 
after reaching a maximum separation move again towards each other. 
This peculiarity seems to be characteristic of the Tradescantia cells, 
so far as Barber’s data go. It may be due to the presence of relatively 
rigid cell walls, which make the dynamics of the cell division different 
from that of a freely elongating cell. On animal cells there are four 
sets of data: two on the first meiosis in the spermatocites of Stenobo- 
thrus, for two different chromosomes in the same cell, and two on the 
second meiosis, for two sister cells. In the first two cases the situation 
is again complicated by the presence of chiasmata. Thus only the last 
two sets seem to fulfill the simple conditions assumed above. 

As seen from Figure 1, both sets of data are represented within 
the limits of experimental error by the simple equation (2). There 
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FIGURE 1 
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seems thus to be no indication in this case of the presence of any re- 
pelling forces. 

This immediately raises the question on the nature of mechanical 
stability of the equatorial plate. Quite regardless of the nature and 
mechanism of the forces present, the chromosomes will be located 
stably in the equatorial plane only under one of the two conditions: 
either they are attracted to the poles by forces which increase with 
distance, or else they are repelled by the poles by forces which de- 
crease with distance. 

The only known forces, that increase with distance are the elastic 
forces. Thus we seem to be led to the conclusion that the elastic spindle 
fibers are not only responsible for the movement of the chromosomes 
at anaphase, but also for their arrangement in the equatorial plane at 
metaphase, unless we assume that some forces are present during 
prophase and metaphase, which disappear at the beginning of ana- 
phase. 

In any case these considerations support quantitatively the theory 
that the movement of chromosomes during anaphase is due to an 
elastic pull by contracting spindle fibers. 


LITERATURE 


’ Barber, H. N. 1939. “The Rate of Movement of Chromosomes on the Spindle.” 


Se we Re, 


Chromosoma, 1, 33-50. 
Rashevsky, N. 1938. Mathematical Biophysics. Physicomathematical Foundations 
of Biology. Chicago: The University of Chicago Press. 
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ON REINFORCEMENT AND INTERFERENCE 
BETWEEN STIMULI 


GALE YOUNG 


THE UNIVERSITY OF CHICAGO 
(AT PRESENT: OLIVET COLLEGE, OLIVET, MICH.) 


eit is shown that the current “two-factor” theory of nerve excita- 
tion can account for sustained inhibition or enhancement by a se- 
quence of stimulus pulses, and for the decrease in the reinforcement 
period with each successive pulse of the train. 


1. Excitation equations. The basic equations in Rashevsky’s exci- 
tation theory are* 


de 

Te hme ae ’ 

dj —Mi=m a 
dt A} s 


which describe the variation of a quantity «=e — 7 under the influ- 
ence of a stimulating current J(t) reckoned as positive at the cathode 
and negative at the anode of the stimulating electrodes. The equations 
are considered to hold so long as o remains below a positive threshold 
value h. When « becomes equal to h, excitation occurs, and the theory 
makes no claim to hold during the succeeding cycle of membrane ac- 
tivity. It is valid again by the time the nerve has returned to its rest-» 
ing state, and presumably somewhat earlier than this, such as near 
the beginning of the relative refractory period. 
The solution of the equations (1) is given by 


Ch er eK [1@erae] ; 
f (2) 
j(t) =e™[jo + M [ I(x) em? da]; 


*Here « and j take the place of ¢ — & and j — j, in one of Rashevsky’s nota- 
tions (1940, p. 108), and are not restricted to positive values. Our threshold h 
corresponds to his quantity j, — ®- The present notation is that used by him in 


another place. (1940, p. 123). 
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where the zero subscripts denote values at t = 0.7 

Equivalent formulations of the theory have been given by Hill 
and by Monnier, but for our present purposes the Rashevsky form 
(1) is the most convenient. 

For application to actual nerves it is appropriate to restrict at- 
tention to the case 


Ko Ma, 
(3) 
ee 0 Ue. 
and without loss of generality we may set 
(mt Me (4) 


2. No stimulus. Starting from an initial state of the nerve 
specified by values e, and 7, it is seen that in the absence of stimula- 
tion, o varies according to 


o = & e*t —j,em , (5) 


Depending upon the values «, and 7, there are various possibilities for 
the time course of «. In any case both o and do/dt vanish at t = o ; 
in the various cases they each have 0 or 1 finite positive roots. The 
zero value of « comes at the time 


1 &5 
t= ie 


provided that this quantity is real and positive; while a maximum or 
a minimum occurs at 


ke, 


k-—™m™ log Mio ” (7) 


tin — 


if this is real and positive. Since k > m it is seen that t,, exists if t, 
does, and that t, > t,. There are three cases as illustrated in Figure 
1, and occurring as follows: 


(a) and t,, exist. 

(b) Le - = - ; tm exists, but not ty . (8) 
0 

(c) < a ; neither ¢,, nor t, exist. 3 


inenecute of the diagrams shows that the only cases in what 


+ Rashevsky’ s €&, and j, (1940, p. 108) have a different meaning. 
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o rises to a positive value greater than o are in (a) with o, negative 
and in (b) with o positive. Upon reference to (8) it is found that 
for one or the other of these cases to exist, it is necessary and sufficient 
that the three inequalities 


e<0 
7<0 (9) 
ke— mj <0 


shall all be initially satisfied. Hence excitation cannot occur in the ab- 
sence of a stimulus unless previous stimulation has brought the nerve 
to a state in which (9) holds. Condition (9) is necessary but not 
sufficient for spontaneous excitation to occur. An example of excita- 
tion occurring in the absence of a stimulus is at the anode after a 
sufficiently strong current has been broken. 

If a nerve free of stimulation does not satisfy (9) at one moment, 
then it never will at any later moment. On the other hand a free 
nerve can satisfy (9) only up to the time ¢,,, since the quantity on 
the left of the last inequality is simply —do/dt and so changes sign at 
e., 


FIGURE 1 


From (2) it is seen that a positive stimulus (J > 0) applied toa 
resting nerve (& —j.—0) can never bring the nerve into a state in 
which (9) holds, and thus excitation cannot occur after the cessation 
of a positive stimulus applied to a resting nerve. 

3. Instantaneous pulse. For t very small the exponential factors 
in (2) are sensibly equal to unity and thus after a very brief or in- 
stantaneous pulse we have 


e=e,+KQ 
j=Jo + MQ ; 


~ where Q= [rea ae (11) 


(10) 


8 MATHEMATICAL BIOPHYSICS 


is the total quantity of electricity passed during the pulse. Upon 
reference to (3) it is seen that a positive or cathodic (i.e. Q > 0) 
pulse increases all three quantities on the left sides of the inequalities 
(9). Hence a nerve which is about to fire spontaneously can be stopped 
from so doing by a sufficiently strong positive pulse. A positive pulse 
cannot leave a nerve in state (9) unless it was already there. 

From (10) it follows at once that 


¢=ar (KW) 0. (12) 


Let Q, denote the positive pulse strength just sufficient to fire the 
resting nerve. Then upon setting o—h —1, o = 0, it is seen that 


ae 
ee (13) 
Qo 
If Q is the positive pulse just sufficient to produce immediate excita- 
tion of the nerve in some other state,* then the value of o in that state 
just before the pulse passed was 


ae. 
Om 

Thus the value of o can be experimentally determined (Hill, 1936; 
Coppee, 1936; Katz and Schmitt, 1940), providing a test for the 
theory. Apart from the constant unity, o is merely the negative of 
the “relative threshold”’ and may thus be called the “excitability”. Its 
range is from —o , corresponding to absolute refractoriness, up to 
unity, at which excitation occurs. 

4. Sequence of instantaneous pulses. After a positive pulse of 
strength Q is applied to a resting nerve, o is given by (5) with 


& = KQ 
jo= MQ. 


o=—1 (14) 


(15) 


Since K > M, we have case (a) in (8) and in Figure 1 (cf. Hill, 1936, 
Figure 76). A second positive pulse following the first within an 
interval up to t, will be reinforced, while if it comes later than t, , it 
will be opposed by the stimulus trace or residue of the first. For any 
two sub-liminal pulses such that Q, + Q. > Q, there will be a period, 
greater than zero and less than ¢, within which they will “sum” to 
bring about excitation. This interval depends on the nerve and on 


* The word immediate is needed to handle those cases in which the nerve was 


about to fire spontaneously. A zero pulse would then be followed by excitation, 
but not immediately. 


ig lll ae ie Dea 
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the strengths of the stimuli. Furthermore, if the strengths are un- 
equal it depends on the order in which they are presented, being great- 
er if the weaker stimulus is presented first. 

At atime A after the first pulse, we have 


e—KOe-* 
j=MQe™ ; a 
and if another equal pulse be then applied these values jump to 
e— KO (1 7 €**) (17) 
j=—=MQ(1+e"). (18) 


Continuing in this manner with equal and equally spaced stimuli we 
have just after the n-th pulse 


e=KQ(1 + e* + e78 $ oF emaiay , (19) 


with a similar expression for 7. Upon summing the geometric pro- 
gressions we obtain 


Bake — 
1 == e-nmd (20) 
j=MQ>—~ 
which, as 7 becomes large, reduce to 
c= KO oe 
: (21) 


Dag A has pcarer-r tl 


It is seen from (20) that « and 7 are always positive. Hence (9) 
is never satisfied and excitation cannot occur between pulses. Also, 
if o is less than zero just after a pulse it remains so until the next 
pulse, as may be seen by noting that case (a) in (8) is not possible 
with 7) > €& >0. Consider in particular the steady state for which 
(21) holds. o will be less than zero at the beginning of a cycle (just 
after a pulse) if and only if the quantity 


K1i-e™ 


5 eR Rennes 22 
Biralisaie es a 
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is less than unity. The second factor has its minimum value m/k at 
A= 0 and increases monotonically to unity as A becomes large. Hence 
if 


(23) 


the quantity o can be maintained less than zero throughout the steady 
state cycle if Ais small enough. 

Conversely, if is greater than zero at the end of a cycle 1t must 
have been so throughout. In the steady state the values at the end of 
a cycle (just before a pulse) are 


ere iL 
Ke =e ee 
24 
Hs: MQ e-m& ies MQ 1 ( ) 
d Nea ead eo eee 


Hence o will be greater than zero throughout the cycle if and only if 
the quantity 
Ne een | 


M1 ee 


is greater than unity. The second factor has its maximum value m/k 
at A = 0 and decreases monotonically to zero as A becomes large. 
Hence if 

ee A 

= Ss = (26) 
the quantity « can be maintained greater than zero throughout the 
steady state cycle if A is small enough. 

In either case, (23) or (26), if A is not sufficiently small then 

o > 0 during the first part of the cycle and « < 0 during the remainder 
of the cycle. The mean value of o over a cycle gives a measure of the 
average increase in excitability. Defining 


A 
As =i ao dt 5 (27) 
0 
we find from (5) and (21) that in the steady state 
nit ee TEL See. 


From (1) it is seen that this state is the steady state value of « which 
would be produced by a constant current of strength I = Q/A, and is 
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thus somewhat of the nature of an electrotonus effect. The relation of 
(28) to (23) and (26) is obvious, and sheds additional light upon the 
Significance of a “normal accommodation” (Hill, 1936; Rashevsky, 
1940, p. 111) described by K/k = M/m. 

From (6) and (20) it is seen that the “reinforcement period” 
after the n-th pulse is proportional to the logarithm of 

Digan bees €7ee 
Mi—e™= {ems - ee 

The last factor decreases monotonically to unity as n increases, and so 
the reinforcement period diminishes with each successive pulse. (Lo- 
rente de No, 1938, p. 240). Since K > M the reinforcement period is 
always positive after the first pulse (n=1), and is in fact equal to the 
time at which o reaches its maximum value when a steady current is 
applied to the resting nerve. If, however, the product of the first two 
factors in (29) is less than unity, then the reinforcement period be- 
comes negative after a sufficient number of pulses and the nerve is 
inhibited from then on as long as the stimuli continue, in accordance 
with the discussion of (22). 

Let o(n) be the value of o immediately after the n-th pulse, and 
denote 

é(n) =e (n+1) —o(n). (30) 


If the first n pulses have failed to fire the nerve, then there is a pos- 
sibility of the next pulse being successful only if 6(n) > 0. From 
(19) it is seen that 


DUG) Cee eM ews), (31) 


so that 6 is negative for n sufficiently large. If 6(1) < 0 then so is 6 
for all larger values of n, and o(”) is a monotonically decreasing 
function of n. If 6(1) > 0 (which, upon reference to (6) and (29), 
is the same as saying that A is less than the reinforcement period 
after the first stimulus) then o(m) first rises and then falls as n in- 
creases, attaining its maximum value at one (or possibly both) of the 
two pulses occurring respectively just before and just after the time 
i K 


pan, 32 
pata (382) 


This time is the reinforcement period for a single pulse and also, as 
mentioned above, the time for the attainment of a maximum o under 
the influence of a constant current. The rise and fall of o under a 
steady current is thus quite closely related to the rise and fall of its 
peak values with a steady sequence of pulses. 
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STUDIES IN THE MATHEMATICAL BIOPHYSICS 
OF DISCRIMINATION AND CONDITIONING I 


H. D. LANDAHL 
THE UNIVERSITY OF CHICAGO 


_ _A mechanism with properties of discrimination and conditioning 
is discussed mathematically with reference to special cases in the 
problem of error elimination: elimination of the longer of two 
paths to a goal, elimination of a blind as well as a return alley, and 
Lashley’s jumping problem. For each case equations are derived 
which are qualitatively correct as far as was determined. Several 
qualitative deductions are made and these are substantiated by 
data available. In principle, the theory makes it possible to predict, 
for any trial, the number of errors at any junction of a maze pro- 
vided certain experimental conditions are satisfied, and if a sufficient 
number of experimental values are given to determine the para- 
meters of the system. 


A theory of error elimination based on a neural mechanism has 
been discussed in previous papers (Rashevsky, 1936, Landahl, 1938b). 
We shall continue the theoretical development by discussing results 
from modified assumptions about the mechanism and its functions, 
and by giving approximate solutions to several specific experimental 
situations. 

1. General theoretical formulation. The basis for the following 
development is the neural mechanism of Figure 1, in which the fibers 
have the properties postulated by N. Rashevsky (1938, chap. xxii). 
The stimulus S, presented alone produces the response R, which, be- 
cause of the nature of the experiment necessarily results in another 
response R,. For example, a rat responds to the sensory cues of an 
alley entrance by entering the alley, and this act leads to some further 
response such as that of eating. The response FR, and the stimulus S, 
may then become conditioned with successive repetitions. Similarly, 
if S,, is presented alone, a response R,,, different from R,, is pro- 
duced, and R,, in turn results in some response RP, , which may or may 
not be the same as R,. S,» may become conditioned to R.. 

When both stimuli S, and S,, are presented simultaneously, then 
under certain conditions a choice must be made, as only one of the 
responses, R, or R,, can be made. When they have both been pre- 
sented, the conditioning centers M and M’ (Rashevsky, 1938, chap. 
xxiv) maintain an excitation acting on C and C’ respectively after S, 
and S,, are discontinued. We assume that C (and C’) represents many 
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FIGURE 1 


circuits having a distribution in thresholds, but that none of them be- 
comes excited unless acted upon by M and V (M’ and VI) simultane- 
ously. This part of the mechanism is the essential part of the mech- 
anism discussed by Rashevsky, (1938) and determines the strength 
of the conditioning [equation (1)] due to additional elements of C 
coming into play with each response R,. The fiber from R,, V, is so 
drawn to indicate that the response R, is in itself a stimulus. 

The centers M and M’ represent many elements having a distribu- 
tion of thresholds. As is implied above, we assume that certain of 
these elements become excited by S, alone, and that number increases 
with S.. But if the state of excitation could be permanently main- 
tained at a constant level, the time, ¢ , between presentation of S, and 
the response R, would not enter into the expression for the strength 
of the conditioned response as it must. Let us therefore assume that 
the centers come under inhibitory influences which thus tend to re- 
move the centers from the steady state of excitation, that is, inhibi- 
tory fibers from various centers act upon M . If these inhibitory ef- 
fects at one of the elements M; exceeds some value h; , the threshold 
of stability of M; , then M; no longer contributes to the excitation act- 
ing on C. Thus the strength of the excitation acting on C decreases 
with the time ¢. Let x be proportional to the random influence (in- 
hibitory minus excitatory factor) acting to disrupt M. If q(x) is the 
probability for a certain x , the chance that those thresholds less than 
h will be disrupted is given by z, the integral of q(x) from h to o. 
But if the value of « changes, on the average, every 7 seconds, then 
the time ¢ required to eliminate those circuits of thresholds less than 
his given by t=7/z. Hence, as z is a function of h, h is a function 
of ¢. But the number of elements acting was originally equal to the 
integral of the distribution function of the thresholds from zero to 
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the maximum (depending on S,). Now after a time, t, the number 
of circuits acting is given by the integral from h(t) to the maximum. 
Thus the number of elements acting is obtained as a function of f, 
given the distribution functions. As the number of elements acting 
determines the strength of the conditioned response, we see that the 
latter is a function of t, [equation (1)] and this relation can be ob- 
tained for any specified case. 

The strength of a positively conditioned response, R, will de- 
pend on the number of trials, n, and upon the time, t, between the 
unconditioned and conditioned stimuli, as well as upon other factors. 
It is almost axiomatic that R will increase monotonically with n, at 
least over a very wide range. N. Rashevsky (1938, chap. xxv) has 
discussed a mechanism of conditioning and obtains R as an exponen- 
tial function of , in which R is zero for nv equal to zero and R in- 
creases to a limit as n becomes infinite. The coefficient of the n in the 
exponent contains the stimulus strength of the unconditioned stimu- 
lus. Assuming that, centrally, this value decays exponentially with 
time, (Rashevsky, 1938, chap. xxvii) we could obtain the function 
R(n, t). More generally we can write for R(n,t) an expression of the 
form 


R(n,t) =¢(t) (L—ewm@). , (1) 


We shall, however, use only the first term of the expansion of 
the exponential in our present discussion. Then, if 


f(t) =¢(%) v(t), (2) 
we have the approximation to equation (1) 
R(n, t) anf). (3) 


Instead of presenting a single unconditioned stimulus, let two 
stimuli be presented, one of which is to be conditioned to the response 
desired, thus referred to as the “correct” stimulus, S., and a second 
stimulus which is not to be conditioned and referred to as the “wrong” 
stimulus, S, (Figure 1). Let ¢ be the cumulative number of correct 
responses or responses to the correct stimulus, and let w be the num- 
ber of wrong responses. For the present, we shall not consider fail- 
ures to respond or other possible responses, so that, if is the total 
number of responses, (hereafter referred to as the number of trials) 


then | 
4 n=ct+w. (4) 


_ The derivative of c with respect to n is the proportion of correct 
responses for the particular values of c and n. Thus, defining P, as 
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the probability of a correct-response for given values of ¢ and n, we 
may identify the probability predicted with the proportion observed 
by writing 


dc 
lee 5 
Pe dn’ (5) 
and similarly 
dw 
— 6 
ae as (6) 


where P,, is the probability of a wrong response. From the nature of 
our restrictions and from an inspection of equations (4), (5) and 
(6), it is evident that P, + P,»=1. 

Let the value of the excitation factor, « , (Rashevsky, 1938, chap. 
xxii) corresponding to the stimuli S, be «.. (synapse s,, Figure 1). 
Similarly let ¢,., be the value of « corresponding to S,. Let A e, be 
the amount of « at synapse s, from center C due to the conditioning 
between S, and the response R,. Also let Ae, be the amount at s, 
from C’ due to conditioning of S,, with R,,. If ¢¢ = &. + Ae, is the to- 
tal amount at synapse s,, and &» = fo + Aéw is the total amount at 
synapse s,,, then, if fibers IJ, III' IV, and IV’, have the properties 
necessary for a discriminating mechanism (Landahl, 1938a), the 
quantity «. — e, plays the role of «, — «, in that development. Using 
those results we have [Landahl 1938a, equation (10)], 


log 2 Py + kee — ey) =0. (7) 


So far we have considered only the possibility of two stimuli be- 
ing responded to. For the more general case of N stimuli, we must 
correct P, for chance by introducing [Landahl, 1939, equation (29)], 


P=n—( P.—= | (8) 


where P is the probability of a correct response corrected for chance. 
Then introducing the above value of P into equations (15) and (10), 


(Landahl, 1938) we obtain instead of equation (7), the more general 
relation 


NP, nea 
Woy | ble exe 08 (9) 


log 


where é,, is the average value of ¢,, for the N—1 wrong stimuli. A dis- 
cussion of equation (9) and a comparison with data obtained for its 
verification will be given in a subsequent paper. For the present we 
shall consider the case for which N = 2, that is, we shall refer to 
equation (7). 
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2. Preliminary applications to maze learning. Let the stimuli 
presented be conditioned to some response which we shall refer to as 
the goal. Let a stimulus S; , an alley entrance, lead to a goal where the 
goal response is elicited. Ri(n;,t:) is the strength of the response 
produced by stimulus S; because of its association with the goal n; 
times. t; is the average time which elapses between presentation of 
S; and the beginning of the goal response. The stimulus S; , after the 
initial trial, is no longer the same stimulus as far as the fibers III, III’, 
IV, IV’ of Figure 1 are concerned, since instead of an initial excitation 
factor «;, being produced, a larger value e; = e€;, + Ae; is produced be- 
cause of the past trials. But Ae; is proportional to R;(n;, t;). Hence 
we may write 


& = & +n Ri (nj, bi) (10) 


where 7 is a proportionality constant. Then introducing equation (3), 
we have 


& =& + Bn f(t), (11) 


where B = ay , n; is the number of responses of S; , and ¢; is the value 
of t when response is made to S;. 

Let a stimulus S, be the entrance to a blind alley. Consider that 
the return response is conditioned to S,, so that the response to enter 
is decreased. (Rashevsky, 1938, chap. xxviii, Landahl, 1938b). The 
values of ¢, , corresponding to S, , will again be conditioned to the goal 
by an amount B 7, f (t,), but at the same time will be decreased by an 
amount F'n, f(t,). Here t, is the time elapsing between presentation 
of S,, and the goal response including the time 2t,,, required to 
travel into the blind and back. The coefficient F corresponds to B , but 
is not necessarily the same, as different responses are being condi- 
tioned. Also, n, is the number of responses to S,. Then we have 


oo Op DN, PUL) F ny f (en) - (12) 


Hereafter we shall assume for simplicity that e.; = &», unless 
otherwise stated. This is equivalent to assuming that responses to S; 
and to S, are equally probable before conditioning has taken place. In 
order to continue the development, it is evidently necessary to use the 
mean values for the various times t. We shall also assume that, on 
the average, the velocity of travel is a constant. This does not at all 
imply that the length of time of running a maze is constant through 
successive trials. 


3. Two routes to a goal response. In this section we confine the 
discussion to the case in which there are two separate paths to the 
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goal. Let a response to one of the paths, designated by the subscript 
c , be called a correct response, so defined because &, < tw, i.e. the cor- 
rect path requires shorter time to traverse. Then from equation (11) 
we have for e, and éy 


ef=%,, + BEIT); (13) 
Ew = Eow 1 Buf (tw) ’ (14) 


where ¢ and w are again the numbers of correct and wrong responses, 
and c + w= 7, the number of trials. The possibility of retracing is 
not considered here. 

Substituting equations (13) and (14) into (7) and recalling that 
Eoc = Eow, and ¢ + w=—7, then introducing (6), we obtain a differen- 
tial equation which may be integrated and solved for w: 


cag 1 
~ BBEF(t.) +f (tw) ] 


W 
2 f(t.) a 


The initial conditions are taken as w = ¢c = n= 0. Unless forced 
trials are given in advance, these conditions are true.. Equation (15) 
is a two parametric curve, the parameters being kBf (t.) and kBf (tw). 
The initial slope is one-half. The curve increases to a limit, which is 
a function of both parameters. In order that the limit exist, we must 
have f(t.) > f(t.). But t. < t», and hence f(t) must decrease with t. 
We shall here consider only the case in which f(t) is given by a 
simple decreasing exponential (Rashevsky, 1938, chap xxviii): 


ACh thee (16) 


where a and g are constants. Then the value of w for n= , w., , May 
be written 


X log 


1 2 e-ate 


Wee By Tay PTA REIT GN LC Lae ere (17) 


The value, w.,, is essentially the total number of errors required in 
learning a task to more or less perfection, an empirical value very 
often used. Expanding (17) we have very roughly, for t,. © t,,, 
[169 =o. Fon te) ] et 
OS (18) 
2kBg 
so that the total number of errors made roughly decreases linearly 
with the time difference in the paths and increases as an exponential 


ee) 


- 
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of the average time length of the paths. Thus, for rapid elimination 
the paths must be as short and as different in length as possible. 

Equation (15) has an upper limit for w. In some instances, at 
least, such will not be the case. If, instead of equation (3) we had 
used equation (1), the resulting w(n) relation would not have had 
this property. The w(n) relation cannot be solved explicitly. For 
large values of w and n, we can obtain the slope (d w/d n),, 


Bk 
dw = ——[o(tc) = O(t0)] 
(aaa 1g) 


Since the slope cannot exceed one-half, we must have b(t.) > $(tw), 
and as t, < t,, then ¢(¢) must decrease with ¢. If this final slope is 
very small, we may use equation (15) throughout with good approxi- 


mation. Gulliksen (1934) has discussed a curve which has an asymp- 
tote which is not horizontal (cf. also Gulliksen and Wolfle, 1938). 


4. Elimination of a blind alley. Consider the case in which two 
entrances are presented, S,, to a blind and S, to the goal. We shall 
not consider retracing so the subject either responds directly to the 
S, and proceeds to the goal, or responds to S,, in which case he goes 
down the blind and back, then necessarily proceeds to the goal. Thus 
retracing must be prevented in an experiment if we are to check the 
following results. Let t, be the time taken in traveling directly to the 
goal and let ¢t, be the time taken to reach the end of the blind. Then 
from equation (11) we have 


Ee = Eve + B(e ne w) f(t) (20) 


as we replace 7 by c and n; by c + w, and ¢; by #, for this case. Whether 
a wrong or correct choice is made, the goal response is eventually ob- 
tained. 

For ¢, we have from equation (12) 


bp eon ly (ty H2kYy = F wf (se) (21) 
as we replace n; by w, subscript 7 by w, and since t;, = t, and 
t; = Ls = Hh . 

Using equations (20) and (21) instead of (13) and (14) we ob- 
tain instead of equation (15) the expression, if «x =F/B, 
1 


Meee B 2a) eT (Gy) 


log {2 f(t) /[2 f(t) — f(4 + 2te) 


+ ef (te) +( FG mdi yast c F(t) Jere Teh 
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for 
fla Hee) ef Gey a208 
or 
1 
was OPEN 9h Ge BGs (23) 
Yaa En ) 
for 


f(t + 2t.) — « f(t) =O. 


Equation (22) is a two parametric equation where the two para- 
meters are k B f(t.) and {k B f(t, + 2t) —k F f(&)}. As the ts 
are experimental values, the only parameter which enters here and 
not in equation (15) is F. If t, in this case is made equal to ¢, in the 
previous equation, then one of the composite parameters is the same 
for the two situations. Whether or not this requirement is so could 
be verified experimentally. 

The limiting value of w, w.., when equation (16) is used, is giv- 
en by 


mat 1 2. e-at1 
Wa Frye g (e-a(Et22) — e-atz) log D e-ati — @-a(ta2tz) 4, e aie 
or (24) 
eats 
W 
2kBg 


corresponding to (22) and (23). In the second equation (24), w. 
evidently increases with ¢t, , but a change in ¢, implies a change in f, 
because of the restriction for this case. In the first equation (24), 
W.. Increases with ¢, if t, is large enough. That is, the number of to- 
tal errors increases for increased distance to the goal, provided that 
this distance is sufficiently large. (cf. Hull, 1934). 

For a constant B , « is proportional to F’,, the coefficient determin- 
ing the rate of conditioning of the return response at the alley end to 
the entrance to the blind. It might be expected that w,, would decrease 
with increasing values of F' and hence x. However, for sufficiently 
large values of t and sufficiently small values of t,, w,, decreases 
with x. 

The variation of w., with t, is complex. For large enough values 
of t,, W.» increases with increased t,. But, if t, is small enough and 
if x is small enough, w,, decreases as t, increases. Thus, if « is suf- 
ficiently small, the elimination of a blind will take more trials for a 
longer blind if the goal is far enough away, but if the goal is close 
enough, elimination of the blind will take fewer trials for a longer 
blind. This result has been mentioned previously (Landahl, 1938b, 
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Figure 2) where a qualitative comparison with experimental data 
(Peterson, 1917) is given. If « has a value small enough for such an 
effect to be present, the effect could be eliminated by increasing x, 
either larger F or smaller B . Decreased hunger would mean decreased 
effect of the goal response and hence smaller B or larger x. Thus less 
hungry animals might well show an opposite effect, eliminating the 
longer blinds with more difficulty even near the goal. Such an effect 
was pointed out (Landahl, 1938b, p. 178) to be a possible contributing 
factor in the data by Tolman, et al. (1930) where such a result was 
found. 

We have been considering only a one unit maze. Each unit of a 
maze can be treated in the same way so that we may obtain the varia- 
tion of w with n at each junction. Two difficulties arise. First, the 
time to the goal decreases with successive trials, and that differently 
for the various junctions in the maze. However, since this time varia- 
tion can be predicted theoretically and determined empirically, it can 
be taken into account. In the second place, there enters another factor 
which may be too large to neglect. An animal which is accustomed 
to mazes may have learned to attribute a positive value to a correct 
path because a blind alley is thus avoided. In this case, conditioning 
could take place not only with the goal response but with the learned 
response due to having successfully avoided a blind. Such an effect 
would be especially noticeable far from the goal and would decrease 
the number of errors from that otherwise predicted. This factor 
would considerably complicate the results. 

Even with these two difficulties overcome, we should predict that 
the last alley of a maze would require the same number of errors for 
elimination regardless of the number of blinds in the maze. If there 
is an inhibitory effect from the conditioning centers C and C’ for one 
junction due to those of other junctions, then an increase in the com- 
plexity of the maze would tend to increase the number of errors for 
elimination at each junction. Data by Peterson (1917) indicate fairly 
strongly that the last blind of a long maze is more difficult to elimi- 
nate than the last blind of a shorter maze. The last two alleys in each 
of two mazes were equidistant from the goal. Using the averages 
from the last two alleys in each case, and using the totals of all ani- 
mals on the two sets of experiments, the following values are ob- 
tained: for the six alley maze, the average number of errors for the 
last two alleys is 24, while for the ten alley maze the corresponding 
value is 32. A rough estimate gives the standard error a value of four, 
so that the critical ratio is about two. 

With the above factors accounted for, we should be in a position 
to predict, for any maze, the number of errors for a given number of 
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trials at each junction of the maze. This would include the number 
of errors for elimination of each of the blinds. In order to attempt 
such a prediction, it would be necessary to know the dimensions of 
the maze and sufficient experimental values to determine the para- 
meters. It would also be necessary to know that such factors as left- 
right bias, directional preferences, etc. be properly eliminated or ac- 
counted for. 

5. Elimination of a return alley. Let us now consider one unit 
of a return maze (Landahl, 1938b, Figure 3). We shall assume that 
experimental arrangements eliminate the possibility of retracing. We 
first restrict our discussion to the situation in which only one of the 
entrances to the return alley can be entered, so that there is but one 
correct and one wrong response possible. Let ¢, be the time taken to 
reach the goal if the correct alley is chosen. Let t, be the time taken 
to return to the start if the return alley is chosen. Then from equa- 
tion (11), we have for «, and «, the expressions 


Ec = Eve + Bie te w) f(t), 
6g= a0 6 whee) 


In each trial, whether the return alley is entered, contributing to w, 
or the correct alley is entered, contributing to c, the correct alley is 
entered and conditioning takes place. Introducing (25) into (7) and 
using (6) then integrating we obtain 


1 2 f(t) 


(25) 


ce Se | : 
OEBFE Th) = Bh) —f tb) FP) en 
(26) 
Using (16) we have for the limit of w, 
ea(titte) 1 x 
pays Ose ee (27) 


Thus the total number of trials for elimination of the return increases 
exponentially with the distance to the goal; whereas it increases with 
increased length of the return alley if this length is long enough, but 
decreases with increased length of the return alley if this length is 
small enough. Thus there is one minimum so that a certain length of 
the return alley will give most efficient elimination. But this length 
may be too small to be detected. This result may seem implausible, but 
it may be pointed out that for t, =0, (no alley) w., = c (i.e., no elimi- 
nation), but as a wrong and correct response cannot be differentiated 
there is no contradiction. In this case we have introduced no new 


symbols not found in equation (15), except the ¢’s which are measur- 
able. 
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If now, instead of preventing entrance to one of the return alley 
entrances, we allow either to be entered as well as the goal entrance, 
but again prevent retracing, we will have the same values for «, and 
Ew aS are given in equation (25). Now, instead of using equation (7) 
we use equation (9) with N = 8, there being three choices possible. 
Assuming that the wrong responses are equally probable so that 


Ew = Ew, we have 


1 


= —________. X 
w—EBIG, + &) 
eo Bf (ts) ca 
pelts eee tint &) 5 of (te tls) rn) 0 
for which, using (16) 
weereleah) 1 3 
We TBG 32am (29) 


Since the same parameters occur in equations (28) and (29) as 
occur in (26) and (27), except for the numbers 2 and 3, we can pre- 
dict completely the behavior from one situation to the other. Whether 
or not the prediction is sufficiently accurate could be tested experi- 
mentally. 

It will be noted that, on the basis of this mechanism, elimination 
is essentially abrupt as compared with the previous mechanism in 
which there was gradual elimination. Here we consider only two pos- 
sibilities, entrance or no entrance. However, it can be seen that the 
probability of a chance fluctuation being able to produce the return 
response before the alley end is reached would be larger with succes- 
sive trials. This point requires a more detailed study. 

6. Application to the Lashley jumping apparatus. Let a pair of 
stimuli be presented such that, if response is made to one after time 
#., a reward is obtained, where as if the other stimulus is responded 
to, after a time ¢,,, an unpleasant response follows. Then we have 


Ec = Ec +B cy (t.), (30) 
Ew — Eow — B" W f (tw) ’ (31) 


where the B’s are in general different from each other and different 
from the B in the previous cases. The sign in equation (29) is nega- 
tive as we assume that the result of the wrong response tends to be 


avoided. 


In this case, for illustration, let us permit eo. and &- to have 
separate values. For simplicity let t, = t,. Setting 
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A —— ek (Eoc - Eow) (32) 


and using (6) and (7) we obtain 


1 2 AB’ ‘ 


words war loqunce: ea: Bie oi ae ai 
DT MAE Bd sng Aaa ae Ba ca 
fOr > 2b , and 
Fe lel Gay 2 (34) 
276.8 yea 


forsbe— dee 

The value of A is half of the reciprocal of the initial slope, and hence 
is readily determined by the initial responses. With A obtained, then 
by equation (32) we have k(éoc — €o.), that is, we have a measure- 
ment of the difference in ¢ values for the stimuli (Landahl, 1939, p. 
163). 

Equation (33) is a three parametric curve, the parameters being 
A, k B' f and k B" f. The value of A is determined from the initial 
slope but the other two parameters are not so easily determined. Al- 
though three parameters give more freedom than we might desire, it 
should be noted that it is not in general possible to fit this curve to 
any set of data such that the initial slope is correct, that the limiting 
value is correct, and that the curve passes through some arbitrary 
point in the middle region of the curve. And even where this is pos- 
sible, there is no assurance that the parameters will be positive. It is 
also necessary that the values of B’ and B” vary appropriately with 
change in reward and punishment. Similarly for the more general 
case in which t, ~ t, , a change in either ¢ , or both t’s should have the 
predicted effect upon the parameters. A complete check of the theory 
requires a large amount of data. So, conversely, if agreement is ob- 
tained, the theory unifies a great number of facts. 

That the general shape of the errors-trials relation is correct may 
be shown by comparison with data. In Figure 2, the experimental 
values are by Gulliksen (1934) for rats number (1) and (13), and 
the curves are graphs of equation (33) with 4d = 1, k B’ f = .0121, 
B" = 0 and A = .63, k B’ f = .0229, B” = 0 respectively. Assuming the 
linear relationships of equation (30) and (31) to be correct, this im- 
plies that the effect of the wrong response on the learning is negli- 
gible compared with that of the correct response. According to Gul- 
liksen’s formulation, the ratios of the effectiveness of the wrong re- 
sponse to that of the correct response are 4 and 1.1 in these cases. It 
should be possible to determine experimentally which predictions are 


~ 
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best. If the present formulation is shown to be wrong in this respect, 
or if it is found not to fit the w(n) curves well enough, the only ap- 
preciable change that can be made is the replacement of equations 
(30) and (31) by more exact expressions. 

We have based the above development upon the mechanism rep- 
resented in Figure 1. However, it is possible to proceed on a more 
formal level (cf. Thurstone, 1930, 1933). A very desirable starting 
point for a development analogous to that given above would be some 
law of the effect of reward and punishment (Gulliksen, 1934) togeth- 
er with the law of comparative judgment (Thurstone, 1927, cf. also 
Gulliksen and Wolfle, 1988). For the case of two choices this can be 
done. But for more than two choices a generalization of the law to 
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the case in which N > 2, (or some approximation) would be desirable, 
and this generalization presents some difficulties. As the slopes of the 
learning curve give proportions of judgments for any trial n, the 
scale values corresponding could be obtained from the generalized 
law, and thus the relation between the scale value and the number of 
trials would be obtained. One should expect that this relation would 
be independent of a large number of variables. If so, application 
could be made to specific cases, as above, obtaining analogous expres- 
sions. The present development is formally more similar to those 
mentioned above than is perhaps apparent. 

A comparison of the above equation (15) with Thurstone’s learn- 
ing curve involving the same variables gives the parameters of the lat- 
ter curve, a, \/m and k the following meanings — a is the reciprocal 
of the initial slope and would be N/(N—1) where N is the number of 
possible paths (two in the case treated) ; m is inverse with the dif- 
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ferences in length of the paths; and k is proportional to the quantity 
gBk. The parameter, «, of equation (16) enters in a rather compli- 
cated manner into this comparison preventing a clear cut separation 
into two parts, but the above is essentially correct. Thus, as it will be 
noticed, m is essentially a parameter of the task while k is a parameter 
of the individual as was defined. To this extent the present formula- 
tion is substantiated by data which bears out such a separation in 
Thurstone’s formulation (cf. Wiley and Wiley, 1937). 

The mechanism of Figure 1 is, of course, too simple, but is still 
able to predict some basic relationships. It would be desirable also 
to have it predict a number of other relationships, particularly the 
properties that are formulated by Gulliksen & Wolfle (1938), such 
that predictions can be made for specific types of stimuli and stimu- 
lus relationships. 

The author is indebted to Dr. H. Gulliksen for reading the manu- 
script, and to Dr. A. S. Householder for reading the manuscript and 
checking the calculations. 

This work was aided-in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 


LITERATURE 


Gulliksen, H. 1934. “A Rational Equation of the Learning Curve Based on Thorn- 
dike’s Law of Effect.” J. Gen. Psy. 11, 395-484. 

Gulliksen, H. and Wolfle, D. L. 1938. “A Theory of Learning and Transfer: I.” 
Psychometrika, 3, 127-149. 

Hull, C. L. 1934. “The Goal Gradient Hypothesis and Maze Learning.” Psychol. 
Rev., 39, 25-43. 

Landahl, H. D. 1938a. “A Contribution to the Mathematical Biophysics of Psycho- 
physical Discrimination.” Psychometrika, 3, 107-125. 

Landahl, H. D. 1938b. “Contribution. to the Mathematical Biophysics of Error 
Elimination.” Psychometrika, 3, 169-180. 

Landahl, H. D. 1939. “A Contribution to the Mathematical Biophysics of Psycho- 
physical Discrimination II.” Bull. Math. Biophysics, 1, 156-176. 

Peterson, J. 1917. “The Effect of Length of Blind on Maze Learning.” Behavior 
Monographs, 3, No. 4. 

Rashevsky, N. 1936. “Mathematical Biophysics of Delayed Reflexes in Connection 
with the Theory of Error Elimination.” Psychometrika, 1, 265-273. 

Rashevsky, N. 1988. Mathematical Biophysics. Chicago: The University of Chi- 
cago Press. 


Thurstone, L. L. 1927. “A Law of Comparative J udgment.” Psychol. Rev. 34, 273- 
286. 
Thurstone, L. L. 1930. “The Learning Function.” J. Gen. Psychol. 3, 469-493. 


Thurstone, L. L. 1988. “The Error Function in Maze Learning.” J. Gen. Psychol. 
9, 288-301. tee 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 8, 1941 


MATHEMATICAL BIOPHYSICS OF CELLULAR 
FORMS AND MOVEMENTS 


ALSTON S. HOUSEHOLDER 
THE UNIVERSITY OF CHICAGO 


Rashevsky’s equations for describing the joint variation of cell 
shape and concentration of a metabolite are discussed. Conditions 
for the existence of non-spherical equilibria and the location of these 
are obtained and involve only two parameters a and A. Sufficient 
(but not necessary) conditions for the stability of these equilibria 
can also be expressed in terms of these parameters alone. Necessary 
conditions involve in some cases a third parameter B. Quasi-peri- 
odie fluctuations about a stable nonspherical equilibrium may occur, 
but only in case B lies on a certain finite range which can be de- 
fined in terms of a and A. 


1. The problem. Rashevsky (1939; 1940, chap. v) has employed 
his approximation method to derive the equations describing the de- 
formation of a cell and the variation of the “average” internal con- 
centration of a metabolite in the particular case where the perme- 
ability of the cell for this metabolite is infinite and where the cell vol- 
ume remains constant. With these equations, he shows that, besides 
sphericity, which is always an equilibrium shape, there will be either 
two others, or no others, according to the relations among the physi- 
cal constants involved. For the special case when the internal diffu- 
sion coefficient is much greater than the external diffusion coefficient, 
more complete information as to the occurrence and nature of these 
equilibria was obtained. 

We wish here to employ a transformation by means of which the 
equations can be simplified and their discussion facilitated. From this 
transformation will be obtained the exact equations, readily solvable, 
for locating these equilibria when they occur, and the exact relations 
which are needed to characterize them. 

The equations in question are the following (Rashevsky, 1940, 
chap. v, equations (2) and (8)): 


de _ 8D; De r8 (De + 2Di) + 2(7:De + 2rDi) re 
ai! re (7,D. + 27,D;) (De + 2D;) 
1 dr,  3RTyu DDC Datel (2) 


7, dt 2My (r.Di+7,D.)(2Di+D.). 2 MT 
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For the notations the reader is referred to the source mentioned. We 
point out, however, that the cell is regarded as a spheroid with polar 
radius 7, and equatorial radius 7.. All other symbols refer to con- 
stants, except, of course, t and ¢c, the latter symbol representing the 
excess of the average internal concentration of the metabolite over 
the external concentration at infinity. 

To simplify these equations somewhat, we introduce the symbols 
pe andy: 


p> = 1112” ’ == 11/T2 ’ (3) 


p being the radius of the sphere of equal volume, supposed to be con- 
stant, y being the ratio of polar to equatorial radius. The introduc- 
tion of x , the cube root of this ratio, is for the purpose of rationaliz- 
ing the differential equations. From equations (3) it follows that 


ea Gh eee (4) 


For further simplification of notation we introduce the following 
symbols 


pie ee ee aed; ep ri) e 
D” 4M De® any (Goa Ts: un Ba 
When these substitutions are made, we obtain the equations 
dx qd 20Avee i 
ee AE eA ee eae TR ayes Sa 
dt wells oa = Edit (6) 
de q Zee Ale tee b 
ea SY 2 Pe ee re se 
Ci 3h (3 (x? +4) ), eh) 


equivalent to the pair (1) and (2). 


2. Determination of the equilibria. An equilibrium configura- 
tion is determined by a pair of values, x and c, or y and c, which 
cause the left members of (6) and (7) to vanish.* It is sufficient to 
know the value of x or y at any equilibrium in order to be able to cal- 
culate the value of ¢, so that we shall speak, for brevity, of the equi- 
librium ratios y , meaning by this the value of the ratio y of the polar 
to the equatorial radius of the cell when the cell is in a state of equi- 
librium. Or we shall speak of a spherical equilibrium, of a prolate 
spheroidal equilibrium, or an oblate spheroidal equilibrium according 
to the geometrical form of the cell in the equilibrium in question. 


* The mathematical theory of equilibrium or “critical” points of 
. . 3 . . . . a s st 
of ordinary differential equations, which is presupposed eee can be found Hi 
various places, among which may be mentioned Bendixson (1901). 
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Since for « = 1, the right member of (6) vanishes independently 
of the value of ¢, it follows that the sphere is always an equilibrium 
shape. The corresponding value of c is obtained by setting the right 
member of (7) equal to zero, with « = 1, and solving for c. This 
gives c = I’. In order to study the equations in the neighborhood of 
this equilibrium we set 


x=1+€é, ec=Irt+y (8) 


and expand in powers of & and 7. We retain only the first order 
terms, and find 
dé q 2aA —3(a+1) a 
di abies (ast Ten aindaoo? 
thes sis AE a 
di Lait)? 


(9) 


[a(a+3)Vé— (a+1)2y] +---. (10) 


The characteristic roots of the approximating linear equations are 


@020A — 3 (ai3+-1) qd 


Binns (Guim Baga oe) 


i 
The only parameter in the original equations (1) and (2) which can 
be negative is gq, and since A and I both contain gq as a factor, their 
signs will be the same as that of qg. Since q > 0 during production, 
q < 0 during consumption of the substance being metabolized, and 
since the conditions for an equilibrium to be stable is that 1, and 4, 
shall both be negative, it is at once evident that when the essential 
metabolic activity of the cell is that of the consumption of some sub- 
stance, the spherical equilibrium is always stable; when it is that of 
production of some substance, then the spherical equilibrium is stable 
whenever 


ex Bh dy, 


2a 2) 


A 


In no case is the occurrence of quasi-periodic oscillations (approach 
to the equilibrium configuration by a spiral path in the (z, eis 
plane) possible in the neighborhood of the spherical equilibrium. Here 
and hereafter we neglect all consideration of limiting cases where 
some relations of inequality among the parameters are replaced by 
equalities. 

The non-spherical equilibria, when they exist, are given by set- 
ting the right members of (6) and (7) equal to zero, with x # 1. 
This gives 
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SLAG + a)e 31)? 3) (13) 
Ce 5 Ca a * Oo 
24° + Qax? +a4+1 2aA x 


By equating the right members of these two equations, and suppres- 
sing all factors which are essentially non-null, we obtain the equation 


b(y) = 27 = 24(A — 1) gt ale 0 (14) 


for the determination of the non-spherical equilibrium ratios. 

When q < 0, then A < 0, and equation (14) has no positive roots. 
Hence, when the essential metabolic activity of the cell is that of con- 
sumption of a metabolite, there can exist no equilibrium state apart 
from the spherical one. There is therefore nothing more to be said 
of the case of consumption, and we now confine our attention to the 
case g > 0, when all parameters are positive. 

Either two non-spherical equilibria may exist, or else there are 
none. Evidently if there are two, the roots of (14) must be positive, 
and therefore A > 1. This condition is necessary, but the necessary 
and sufficient condition for the roots of (14) to be real and positive is 
this and the added condition that the discriminant be positive. It fol- 
lows that two non-spherical equilibria exist when and only when 


a+ V2(a +1) 
: : 


A> 


Otherwise there wre none. 

If ¢(1) < 0, then two non-spherical equilibria exist, one being 
oblate spheroidal, and one prolate spheroidal. But this occurs when 
and only when (12) fails. Hence we can say that the spherical equi- 
librium is always stable unless there exist both a prolate spheroidal 
and an oblate spheroidal equilibrium. If (12) is satisfied and there 
exist two non-spherical equilibria, then both are oblate, or else both 
are prolate. In the former case, both equilibrium ratios are exceeded 
by unity, and hence their product is less than unity. In the contrary 
case, the product exceeds unity. From (14) it appears that these two 
conditions are equivalent to a < 1 and to a > 1, respectively. Hence, 
when two non-spherical equilibria exist, one of these is oblate, and the 
other prolate unless (12) is satisfied. But if (12) is satisfied, then 
both are oblate when a < 1, and both are prolate when a > 1. It is 
important to note that the existence and location of the non-spherical 
equilibria is determined by the physical parameters in the combina- 
tions a and A alone, as these are defined by (5). A priori, when other 
things are held fixed, a high metabolic rate, or a large cell size, or a 
high temperature will always assure the existence of these equilibria. 

The conditions for the occurrence of the various cases is illus- 
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FIGURE 1 
This figure represents qualitatively the regions of the (a, A) plane with 
which the various equilibria are associated. A scale drawing of this size fails to 
separate visibly the curves on a drawing. 


trated in the figure. Two of the curves represented are branches of 
the curves 


4(a +4)? — 2(4A—1)?7=1, og — 2A 8 0 
a=1/a. 
The other will be discussed later. The first curve is a hyperbola whose 
center is at (— 1/2, 1) ; the second is a straight line which is tangent 


to the hyperbola at the point (1, 3). The two non-spherical equilibria 
exist and are of opposite types when the point (a, A) lies,above the 


(15) 
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straight line; they do not exist when the point lies below the upper 
half of the right-hand branch of the hyperbola. When the point lies 
between the straight line and this arc of the hyperbola, the equilibria 
exist and are alike in type; both are prolate spheroids when the point 
lies below the point of tangency, both are oblate when the point lies 
above the point of tangency. 

3. Conditions for stability and periodicity. Though the para- 
meters a and A alone determine the location of the non-spherical eqi- 
libria, they do not in all cases determine whether they are stable or 
unstable. To investigate the character of these equilibria, we shall 
expand the equations about each of the equilibrium values. It is in- 
convenient, however, to express these equilibrium values in terms of 
aand A. Consequently, supposing the equilibria to exist, we express 
A in terms of a and the value x, of x at one of them. If 2, is the 
value of x at one of these equilibria, then y, = w,° satisfies (14), and 
therefore 


= 205° + 2002 + a+ 1 


= (16) 


2aA 


Since the product of the roots of (14) is (a + 1)/2, therefore the 
value of x, at the other equilibrium satisfies 


2d i Ole A (17) 


From (18) we have, then, that when x = 2,, 


eeee: 3 (ey a). a 2 
* 24.8 + 2an2 +aFi- (18) 


Hence if we set 
=%+E, c=a4+7 (19) 
in equations (6) and (7) and expand, we obtain 


a ant — 1) 
dt 9BI? (a3 + aja, 


X [-89' (a,3 — 2a) E+ (2a,° + 2an,2 +a+1)y]t.-, 
eS a er ee 
dt = 30x, (443 + a) (20,° + 2ax,5 +a +4 1) a Ot 
+ 4(@ — a — 1a, —a(a+1)}& + (2a,° + aa, (21) 


ee et a Sato 


where only terms of first order have been given. The characteristic 
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roots, 4, and 4, , of these equations have the sum and product 
dat As 


ma —q[ Ba, (22,° + 2ax,3 +a+1) + (#3 —1) (43 — 2a) ] 


SRL? Gy a) pA 


7 (a — 1) (22,* —a— 1) 
7k / ss See ae Se ee 
mids 3BI'x,? (a3 + a) : (23) 


and the discriminant of the characteristic equation is 
(A, — do)? = q?{ [Ba, (2018 + 2aa.3 +a+1) + (43-1) 
X (4,3 — 2a)]? — 12Bx, (x12 — 1) (2,3 + a) (22,° — a (24) 
lf ebe! AC Aa)? 


The equilibrium is stable if and only if 24, + 4% <0 and 4,4. > 0, 
while periodic or quasi-periodic oscillations occur in the neighborhood 
in case it is also true that (4, — 4)? < 0. We have therefore to in- 
vestigate the circumstances under which, at any non-spherical equi- 
librium x, , these inequalities hold. 


Define 
pees (a) il) (7,° — 20) 
eae faint onsale (25) 
a ae 12 C (a+ a) (22,° — @=—"1) 
HOVS(C+ I~ Gam oa) Ga team tarty 


Then a little calculation shows that the condition 1, + 4, < 0 is equiv- 
alent to the condition B > B,, while the condition (A, — 4.)2 < 0 is 
equivalent to f(C) < 0. We note that f(C) becomes positively in- 
finite with C, and that it is positive for C = 0. Hence either f(C) 
> 0 for all C, or else there exists a pair of values C, and C., such that 
for C, < C < C,, it is true that f(C) < 0, while for all other values 
of C,f(C) = 0. Hence either (4, — 4.)? > 0 for all B, or else there 
exists a pair of values B, and B, such that for B, < B < B,, it is true 
that (A, — A.)? < 0, whereas for every other B, (A, — d4.)? > 0. We 
note, then, that f(C) has a minimum at 


10x? = 140n,°-- (40 — 7a — Dee — 4a(@ +1) 


28 
(ee) (Ane ca ee) (23) 


C=C,=+ 


and that the minimal value is 
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HVS 12 (2 Fa) (22° — a — 1) 
Fin = (22 — 2a)? (2x, + 2av2 + a+ 1)? 


| (29) 
x [404° + 8ax,° + 4(a? -—a—1)43—a(at+1)]. 


In discussing the signs of these quantities, we shall thus be concerned 
with the signs of the quantities 


dé: (#1) = 102,9 + 14an,° + (4a? — Ta — Tete = AGE IDE, (30) 


and 


bo (ti) = 40," oda, + 4(0° Oo 1) oe), (31) 
and we note for future reference that 
by Clg) = bay) YE Se ONE res O18 (32) 


We now consider the equilibria in their various possible locations. 
4. Case of two oblate spheroids. We have seen that this can 
arise only when 


Cale (33) 
Let the equilibria E, and EH’, be given by x, and x, where 
conte CO < Pepale. 
Hence, from (17), it follows that 
23° <0 Lae G51 < 246° = 27 (34) 


To examine the equilibrium F,, we have only to replace x, by x, in 
the relevant formulae. From (23), therefore, the interchange having 
been made, it is evident that at H., 4, 42 < 0, and the equilibrium is 
unstable. . But at H;, 4, 42 > 0, and hence the necessary and sufficient 
condition for stability is that 2, + 4. < 0, and hence that B > B,. If 
v3 < 2a, then B, <.0 and necessarily B > B,. But if x3 > 2a, then 
B, > 0, and the equilibrium at H, is unstable unless B > By. 

The condition x, < 2a can be expressed in terms of a and A, if 
we refer to the curve 


(a + 1/2)?=4(A — 47/16) , a=1/a, (35) 


which is the condition that (2a) = ¢(2/a) = 0. This is evidently a 
parabola which opens upward, and is the third of the curves of the 
figure. It is tangent to the hyperbola at the point (—1/2 + vV33/2, 
5), and it intersects the straight line at (2, 9/2) and at (3, 6). When 
the point (a, A) lies above the parabola, then y, < 2a < y,. When 
the point lies between the parabola and the hyperbola, we have either 
2a < Yi < ¥2, or else y; <y2 < 2a, according as ¢'(2a) is negative 
or positive. But ¢’(2a) = 2a(5 — A). Hence when A < 5 and the 
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point (a, A) lies between the two curves, then y, < y < 2a, and 
when A > 5 and the point lies between the two curves, then 2a < y; 
< y2. Hence the equilibrium at EF, is necessarily stable when the 
point (a, A) lies above the parabola, and also when the point lies 
below the line A = 5. In other cases, for points above A = 5, but be- 
low the parabola, the equilibrium is stable only when B exceeds the 
positive quantity B,. 

When two non-spherical equilibria exist and in both states the 
cell 1s an oblate spheroid, a situation which can arise only when 
1/a =a <1, then the spherical equilibrium is always stable, the more 
nearly spherical of the two non-spherical equilibria is necessarily un- 
stable, and the other equilibrium is always stable provided the ratio 
x1° of its polar to its equatorial radius is less than 2a. The condition 
xi° < 2a, when expressed in terms of aand A is that the point (a, A) 
shall lie either below A = 5 or else above the parabola (35). If x° 
> 2a, then the equilibrium at x, is stable only when B exceeds a cer- 
tain positive quantity By defined by (25). The parameters a = 1/a, 
A and B are all defined by (5) and the value of y1 = x,3 at the equi- 
librium is the smaller of the two roots of (14). 

We consider next the possible occurrence of periodic or quasi- 
periodic fluctuations about this equilibrium in case it is stable. If 
x3 > 2a, and B has a value close to B, , then evidently (4, — 4.)2 < 0 
and the solution of the differential equations are quasi-periodic (and 
damped) in the neighborhood of the equilibrium. But if 2,3 < 2a, 
then since (4, — d2)? is certainly positive for B very small, it follows 
that values of B will exist such that quasi-periodicities occur only in 
case the numerator in (24) has a negative minimum for a positive 
value of B. This is equivalent to the statement that f(C) in (27) 
has a negative minimum for a positive value of C. But the values of 
C and of f at this minimum are given by (28) and (29), so that C, 
> 0 if and only if ¢:(#1,) < 0 while fnin < 0 if and only if ¢.(x.) < 0. 
But from (382) it follows that the second relation implies the first. 

Consider, therefore, the inequality 


Ay,® + 8ay,? + 4(@ —a—1) Yy, —a(a+1) <.0,; (36) 


where y, = 21°, with the assumption y, < 2a. We have seen that for 
this it is necessary that a > 1/3. For y, = 2a, the left member of 
(36) has the value 9a(8a? — a — 1). Consequently if 8a? -a—1< 
0, (36) is always satisfied for y, < 2a. But if 8a2-a—1> 0, (86) 
_ will be satisfied only by values of y, , less than a certain value y’; < 2a. 
Hence the situation may be described as follows: 

When two oblate spheroidal equilibria exist and y, = x,° repre- 
sents the smaller equilibrium ratio of polar to equatorial radius, then 


36 MATHEMATICAL BIOPHYSICS 


whenever 8a2 — a — 1 < 0, there always exists a pair of numbers B, 
and B. such that when B satisfies B, < B < Bz the equilibrium yi 8s 
stable and in the neighborhood of this equilibrium quasi-periodic fluc- 
tuations in y and in c occur as the cell approaches equilibrium. But 
if 8a2 —a — 1> 0, such values B, and Bz exist only in case y, is small 
enough to satisfy (36). When 8a2 — a — 1 > 0 and (36) fails, then 
whatever value B may have, the equilibrium will be stable and no fluc- 
tuations will occur. 

5. Case of two prolate spheroids. We have seen that this can 
arise only when 


fia (37) 
Let the equilibria LH, and E, be given by x, and 7, where 
bs Ye Se 
Hence it follows that 
Oe Ol ae e Dates <n ee neha (38) 
Since a > 1 it follows that 4a > a + 1, whence 
Gao Ase (39) 


To examine the equilibrium FE, we replace x, by x, in the relevant for- 
mulae. From (23), therefore, it is evident that at E., A. 4. < 0 and 
the equilibrium is unstable. But at EH, , 4, 42 > 0, and hence the neces- 
sary and sufficient condition for stability is that 4, + A. < 0, and hence 
that B > B,. It is evident from (25) that B, > 0, and from (24) 
that for B= By, (dA, — 42)? < 0. Therefore we can say 

When two non-spherical equilibria exist and in both states the 
cell is a prolate spheroid, a situation which can arise only when a > 1, 
then the spherical equilibrium is necessarily stable, the more nearly 
spherical of the two non-spherical equilibria is necessarily unstable, 
and the other equilibrium is stable provided the parameter B exceeds 
a certain positive quantity B,. If B also exceeds a certain quantity 
B. > Bo, then a configuration in the neighborhood of this equilibrium 
approaches the equilibrium monotonically. But if By) < B < B., then 
the cell, when brought into the neighborhood of this equilibrium will 
approach the equilibrium through quasi-periodic oscillations about 
the equilibrium state. Criteria for these circumstances depend eu- 
clusively upon the combinations a, A , and B of the physical constants 
as defined in equations (5) ; the ratio y, = x,° of the polar to the equa- 
torial radius is defined by (14), of which it is the larger root, Bo is 
defined by (25), and B, is the larger of the two values of B which 
cause the right member of (24) to vanish. © 

6. Case of one oblate and one prolate spheroid. Let the equi- 
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libria EL, and E, be given by x, and x, where 


io Pe (40) 
so that, by this and (14) it follows that 
20° < 24> << 4+ 1, Grae <n (41) 


For discussing the equilibrium of the prolate spheroid, x, must be re- 
placed by x2 in the relevant formulae. From (23) and (41) it is evi- 
dent that 4, 4. > 0 at both equilibria. Then (25) defines Bo; corre- 
sponding to EL, and corresponding to H, we may define 


PRE 1) la 20) 


i. = eee eee 
0 ea (2x—® + 2ax2 + a+ 1) 


(42) 
The equilibrium E, is stable provided B > B,, and E, is stable pro- 
vided B > B,.. We note that B,, > 0 when x,’ > 2a, and Bx. > 0 when 
v2 < 2a. Hence B,, and B,2. cannot both be positive because of (40). 
From (14) we see that B,, and B,. are both negative and both equilib- 
ria are necessarily stable when the point (a, A) lies above the curve 
(35). If a < 2 and (a, A) lies below (35), Ba < 0, Boo > 0, so that 
while the oblate spheroid is necessarily stable, the prolate spheroid 
may be unstable. If a > 3 and (a, A) lies below (35), the prolate 
spheroid is necessarily stable but the oblate spheroid may be unstable. 
Let there exist two non-spherical equilibria, one in which the 
cell is an oblate spheroid with the ratio y; = x,', of polar to equatorial 
radu, and one in which the cell is a prolate spheroid with the ratio 
Yo — x.*. These ratios are determined by the equation (14), in which 
the a and A are defined by (5). The spherical equilibrium is always 
unstable, whereas at least one of the non-spherical equilibria is bound 
_ to be stable: the equilibrium x, if the point (a, A) les in the lower 
region between the curve (35) and the second of curves (15), the 
equilibrium x, if the point lies in the upper of these regions, and both 
if the point lies above both curves. In either of the first two cases, the 
_ other equilibrium is also stable provided B exceeds the positive one 
of the two quantities B,. as defined by (25) or Boz as defined by (42). 
j Consider finally the possibility of quasi-periodicities. If B,. is 
positive, the cell will approach the equilibrium state x, through 
quasi-periodic fluctuations, provided B is sufficiently close to Bo. If 
Bi. is positive there will be a similar approach to x, provided B is 
sufficiently close to By» . 
Suppose B,, < 0. Periodicities exist for suitable values of B pro- 
_ vided f(C) as defined by (27) has a negative minimum for a positive 
- value of C. But B,, < 0 only when x,° < 2a, and hence for periodic- 
ities we must have ¢,(%7) < 0 and ¢.(71) < 0, where ¢: and ¢2 are 
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defined by (30) and (31) respectively. But, in view of (32) it appears 
that $2(2:) <0 implies ¢,(a:) <0. If 8a? — a — 1 < 0, then $2 (#3) 
<0 for any 2,° < 2a, since ¢2 < 0 when #,° = 2a. If 80° — a — 140; 
then a fortiori a > 1/3, whence 4a > a + 1h 2058 - When 2,3 = 
(a + 1)/2, then ¢2(%1) =38(a + 1) (8a? — 1)/2. Hence if 8a? — Gil 
> 0 but a < V3/3, ¢o(a1) < 0 necessarily. But if a > V 3/3 then 
there is a single positive a'13 << (a+ 1)/2 < 2a such that ¢2(2'1) = 90, 
and if 7, < #1, ¢2(%) < 0. 

In the case of a single oblate spheroid of equilibrium with x,° the 
ratio of polar to equatorial radius, there is at most a finite range of 
values such that if B , as defined by (5), is on this range, there will be 
quasi-periodic fluctuations in the vicinity of the stable equilibrium. 
This range always exists unless a > V3/3, which implies that x° < 
2a, and it will exist also in the latter case unless x, is such that 
2(Xi) > 0, where d. is defined by (32). 

At the equilibrium 2, , if Boo > 0 a range of values of B always 
exists for which quasi-periodicities exist. We can have B,. < 0 only 
when «7,2 > 2a, as appears from (42). Then the conditions for the 
existence of such a range are ¢:(%) > 0 and ¢2(%,) > 0. But from 
(32) and (41) it appears that the second of these inequalities implies 
the first. If 8a? — a — 1 > 0 then the second inequality is satisfied 
for any x,° > 2a since it is satisfied when x, = 2a. If 8a2 —-a — 1 
< 0, then there will exist an wx’, such that w’,° > 2a, d2(x’2) = 0, and 
when 2, > x’. then such a range exists. 

In the case of a single prolate spheroid of equilibrium with xo* 
the ratio of polar to equatorial radius, there is at most a finite range 
of values such that if B , as defined by (5), is on this range, there will 
be quasi-pertodic fluctuations about the stable equilibrium. This range 
always exists unless 8a? — a — 1 < 0 and x,’ > 2a, and if this is true 
the additional condition $2(X2) > 0 must be satisfied in order for the 
range to exist. 

The author is indebted to Dr. A. Weinberg for reading the manu- 
script and for his valuable suggestions. 
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